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Abstract 

We examine optical F-band light curves in luminous eclipsing black hole X-ray binaries, using a super- 
critical accretion/outflow model that is more realistic than the formerly used ones. In order to compute 
the theoretical light curve in the binary system, we did not only apply the global analytic solution of the 
disk, but also include the effect of the optically thick outflow. We found that the depth of eclipse of the 
companion star by the disk changes dramatically when including the effect of the outflow. Due to the effect 
of outflow, we could reproduce the optical light curve for typical binary parameters in SS433. Our model 
with an outflow velocity v ~ 3000km/s could fit whole shape of the averaged F-band light curve in SS433, 
but we found a possible parameter range consistent with observations, such as M ~ 5000-10000 L^/c 2 
(with Le being the Eddington luminosity and c being the speed of light) and Tc=10000K-14000 K for 
the accretion rate and donor star temperature, respectively. Furthermore, we briefly discuss observational 
implications for ultr aluminous X-ray sources. 
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1. Introduction 

The profile of a light curve in an astronomical object 
includes a wide range of information, i.e., the spatial dis- 
tribution and dynamical motion of gas, or its radiative 
processes, etc. The analysis of the light curve has been 
therefore known as an important tool from the old days. 
Especifically, in the case of the eclipsing binary system, 
which has a compact star, it is useful to suppose the 
brightness distribution of the accretion disk around the 
compact star. This method can also be applied to black 
hole candidates. Generally, it is difficult to discern light 
from the companion star and light from the accretion disk, 
but if the object has an eclipse, it may be possible for 
two ingredients to separate. Moreover, black hole candi- 
dates show the (partially) eclipsing property, and we could 
close in on the feature of the accreting gas, or the black 
hole mass, or spin, etc (Fukue 1987; Watarai et al. 2005; 
Takahashi & Watarai 2007). 

For the last decade, a large number of luminous 
black hole candidates (BHCs) have been found in nearby 
galaxies. They are called "ultr aluminous X-ray sources 
(ULXs)", but their origin is still unknown (Makishima et 
al. 2000; Roberts et al. 2005). X-ray data analysis is 
the mainstream study of the ULXs, but other band data 
will be useful for evaluating other binary parameters. It is 
no wonder that eclipsing binaries are detected among the 
BHCs. In fact, several eclipsing black hole binaries have 
been discovered (Orosz et al. 2007 in M33; Ghosh et al. 
2006 in NGC4214). 

To examine the observational properties of such lu- 



minous eclipsing binaries, light-curve fitting will be a 
powerful tool. The conventional study of a light curve 
in a binary system applies a geometrically thin disk, a 
so-called standard accretion disk (Shakura & Sunyaev 
1973). This model may have applicability to the high/soft 
state in black hole candidates, but it may not extend 
to the low/hard state, very high state, or super-critical 
state, which is a more luminous state than the high soft 
state. Luminous black hole candidates that exceed to the 
Eddington luminosity seem to accrete a large amount of 
gas via a companion star, and the disk becomes geomet- 
rically thick. This type of accretion flow is called a "slim 
disk" (Abramowicz et al. 1988; Watarai et al. 2000). The 
geometrical thickness of this disk can play an important 
role in covering a companion star, depending on an angle, 
and contribute to changing the shape of the light curve. It 
is necessary to include it in a calculation of bright binary 
system properly. Fukue et al. (1997, 1998) performed 
light curve analysis in SS433 using geometrically thick 
torus model by Madau (1988), but it is known that the 
disk model is thermally unstable. Hirai & Fukue (2001) 
compared the optical light curve in SS433 with (thermally 
stable) supercritical accretion disks, but they adopted the 
self-similar solutions even at the disk outer region. We 
thus adopted a thermally stable, more realistic (appropri- 
ate) treatment of the disk model, and compared the model 
with observations of super-critical accreting objects such 
as SS433 and ultr aluminous X-ray sources. 

In addition, in the previous studies the mass outflow 
from a supercritical disk was not included, but a naked 
supercritical disk was considered. Hence, in the present 
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study we consider the supercritical accretion disk with an 
optically thick outflow from the center. 

In the next section, we introduce the assumptions of 
our model. In section 3, we briefly show the light curve 
calculation method. Calculation results are presented in 
section 4. The final section is devoted to concluding re- 
marks. 

2. Model for Accretion and Outflow under 
Supercritical Accretion Flows 

In this section, we calculate the light curve when the 
binary star system that fills Roche-lobe is assumed, and 
a supercritical accretion has occured in the compact star. 

Ideas of supercritical accretion have been proposed 
by many authors in early 1970's (Shakura & Sunyaev 
1973; Abramowicz et al. 1978; Jaroczynski et al. 1980; 
Paczyhski & Wiita 1980). Whether accretion that ex- 
ceeds the Eddington luminosity is possible or not has 
been doubted for many years. The solution for supercriti- 
cal accretion had already been obtained in one dimension 
(Abramowicz et al. 1988). Due to the development of 
the computer, supercritical accretion is actually known 
to be reproduced by numerical simulations (Okuda 2002; 
Ohsuga et al. 2005; Ohsuga & Mineshige 2008). It thus 
becomes impossible for us to deny supercritical accretion 
any longer. In the next subsections, we briefly introduce 
our model and assumptions. 

2.1. Model for Supercritical Accretion Flows 

The photon trapping radius characterizes supercritical 
accretion at the radius where advective energy transport 
becomes important (Begelman & Meier 1982; Ohsuga et 
al. 2002). This radius accords with the radius that the 
gravity of the disk balances with the radiation pressure, 
and it has been suggested that an outflow may originate 
inside it (Shakura & Sunyaev 1973; Lipunova 1999; Fukue 
2004; Heinzeller & Duschl 2007). 

However, the outflow blows a very tiny area compared 
with the size of the whole disk. It actually will not have an 
influence on computing optical flux. Recently Takeuchi et 
al. (2009) analyzed 2D RHD simulation results by Ohsuga 
et al. (2005), and rebuilt a one dimension model, but 
the effective temperature distribution of the one dimen- 
sion model hardly changes. Therefore, we can use an one- 
dimension model safely without any outflow effect. 

2.1.1. Radiation-pressure dominated regime: ft es ^> 

For radiation-pressure dominated regime, Watarai 
(2006) constructed analytical formulae that can be ap- 
plied for a wide range of accretion rates, and it has been 
shown that these solutions are a good approximation of 
the numerical solutions. According to Watarai (2006), the 
scale height of the disk is given by 

H & = 3.0f(r,m) 1 ' 2 r. (1) 

This solution is characterized by the ratio of the advec- 
tive cooling rate to the viscous heating rate, i.e., / = 
Qadv/^vis' wmcn can be represented by an analytical 



form dependent on the radius and the mass accretion rate, 
/(r,m). The radius f is normalized by the Schwarzschild 
radius, and the m represents the mass accretion rate in 
Eddington units (M~Edd = L^/c 2 ). The explicit form of 
/(r,m) is given by 

1 



D - 2(f/m) 2 + 2 - D(r/m)y/D 2 (r/rh) 2 + 4 (2) 



where D is a constant of order unity (e.g., D « 2.18 for a 
polytropic index N = 3). The function /(r,m) is close to 
unity for an advection dominated regime, and it is close to 
zero for a radiative cooling dominated regime (see Watarai 
2006 for more details). 

The effective temperature distribution is given by 



T eff ^4.48xl0 7 /(r,m) 



1/8 



m 



-1/4^- 



X / 2 K. 



(3) 



The boundary radius between the radiation-pressure 
dominated regime and the gas-pressure dominated regime 
is located at 

rrad-gas = 18 (am) 2 ' ' 3 m 16 ' ' 21 (4) 

* 601(a/0.1) 2/3 (m/10) 2/3 (m/100) 16/21 . (5) 

The same equations are posed by Shakura & Sunyaev 
(1973). We note that analytic solutions shown here are 
useful for radiation-pressure dominated regime. 

2.1.2. Gas-pressure dominated regime 

In the gas-pressure dominated regime, if electron scat- 
tering dominates the opacity, the scale height of the disk 
and the effective temperature distribution are 



H b = 2.7 x lO'a-^m^m^f 21 / 20 , 
T eff « 3.50 x lOFm-^th^f-^K, 



(6) 
(7) 



with the same formula by Shakura & Sunyaev (1973). 
The transition radius where ft es ~ ^ff is 

r g as,out = 2.5xl0 3 m 2 / 3 , (8) 

and thus it is a simple function of the mass accretion rate. 

The outer region of the disk is dominated by free-free 
opacity, and the scale height is given by 



H c = 1.5 x lO 3 ^- 1 / 1 ^ 9 / 10 ^ 3 / 20 ^/ 8 , 



(9) 



with the effective radial dependence of the temperature as 
in equation (7). 

We ignore the irradiation by the disk itself or the photo- 
sphere of the outflow, because the irradiation dominated 
regime appears at the outer region of the disk, and the 
temperature and geometrical effects do not contribute to 
the optical light curves (less than 10 %). To avoid confu- 
sion of the model, we decided to handle an easier model. 

2.2. Model for Massive Wind 

Shakura and Sunyaev (1973) proposed the massive (su- 
percritical) outflow, which is formed by the strong radi- 
ation from the disk. The size of the photosphere surface 
made by the outflow becomes large and has an influence 
on the form of the light curve. Thus, we should include 
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the effect of the outflow in our model. We do not include 
the collimated, relativistic jet component in this paper. 

Here we introduce a simple wind model by Abramowicz 
et al. (1991), which assumes a spherical symmetry and 
uniform outflow velocity. The density of the wind, p w , is 



M 



47T^7 



R~ 



(10) 



where M oll t is the mass outflow rate, v is the velocity of 
the gas, 7 is the Lorentz factor: (1 — /3 2 ) 1 / 2 , where f3 = v/c, 
and R is the distance from the black hole. 

The mass outflow rate should be determined by the 
physics of the interaction between the disk and outflow, 
i.e., M ut = vM acc , where n is the efficiency of the outflow 
gas from accreting gas. In our present study, we assume 
r] = 1.0 for simplicity. That is, our model simply assumes 
that all accreting matter changes to the outflow at the disk 
inner edge (3r g ). The location of the boundary layer be- 
tween the outflow and the disk may have a strong impact 
on the emerging X-ray spectrum, but it is not expected 
to have an enormous influence on the optical band, since 
the optical emitting region of the photosphere is far away 
from the disk inner edge. 

3. Binary Light Curve Calculation 

To calculate the F-band flux in a binary system, we 
adopt the "Ray- Tracing Method". We suppose that the 
binary star fills the Roche lobe and transports its mass 
to the compact star. The shape of the companion star 
reflects the shape of the potential. The photon propa- 
gates from an emitting point on the surface of the disk 
or that of the photosphere of the wind to the distant 
observer. According to Fermat's principle, however, the 
light rays are traced from the observer's display coordi- 
nate (x s ,y s ,z s ). After the ray arrives at the surface of the 
disk/wind, we evaluate the geometrical thickness and ef- 
fective temperature of the disk/wind model presented in 
the previous section. The observed flux is integrated in 
each optical bands (V, R, I), assuming the blackbody ra- 
diation at the surface of the disk, star, and photosphere 
of the outflow. 

The spatial resolution of the calculation is about 1% — 
0.1% of the binary distance, which equals to 10 4 r g - 10 3 r g . 
This resolution is sufficient to resolve optical flux from the 
binary system. Thus, mesh size of this calculation is at 
least much smaller than 10 % of this radius. 

3.1. Location of the Photosphere 

Location of the wind photosphere is estimated from the 
point where the optical depth measured from a distant 
observer equals to unity, that is, 
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7(1 — f3 cos 0)hvp w ds. 



(11) 



ph 



where the i? p h is location of the last scattering surface of 
the photon in the disk coordinates, 7 is the Lorentz factor, 
Kj is the opacity for electron scattering, f3 is the velocity of 



the wind normalized by the speed of light c, and the 6 is 
the inclination angle. Abramowicz et al. (1991) obtained 
an analytic solution of a moving plasma, and the shape 
of the photosphere in their model does not have spherical 
symmetry. We use their formulae in this paper. 

In some cases, analytic formulae derived by King (2003) 
are useful to estimate the size of the photosphere. It is 
given by 

^ ph ~i^r' (12) 

where i? p h is the radius of the photosphere, k is the opac- 
ity for electron scattering, and b is the eccentricity of the 
photosphere. The value of b is 0.5 — »• 1, which is fixed at 
unity throughout this paper. 

In this paper, we assume that the mass outflow rate, 
Moutj is equal to the mass accretion rate at the central 
region. The mass outflow rate in our model gives the 
maximum rate. Typical size of the photosphere is given 

by 

Rph _ T^out _ -^q5 ( ^out X 



2/3 



2000 



P 
0.01 



(13) 



As we show later, if the size of the photosphere is much 
smaller than the disk size, the photosphere does not influ- 
ence the shape of the eclipsing light curves. 

3.2. Temperature and Luminosity at the Photosphere 

We evaluate the maximum temperature of the photo- 
sphere by using the following procedure. First, we esti- 
mate the size of the photosphere of the outflow. Assuming 
that photons ejected from the disk surface inside the pho- 
tosphere is conserved until escaping photons from the sur- 
face of the photosphere, i.e., all photons are generated 
inside the disk, not be generated in the wind. The tem- 
perature of the wind photosphere T w then is given by 



^disk 



47r(6i? ph ) 2 • 

Here the luminosity of the disk, Ldisk, is given by 
W= / 2<rT e 4 ff • 2nrdr. 



(14) 



(15) 



To determine the location of the disk surface or the pho- 
tosphere of the outflow, we integrated the optical depth 
from an observer at infinity to the surface of r p h = 1 along 
the line of sight. After the light rays arrived at the sur- 
face, we measured the temperature using equation (14) 
and the geometrical thickness of the disk from the under- 
lying disk/ wind models. 

4. Optical Flux Images 

4-1. Case without Outflow 

Let us see the results when only a naked supercriti- 
cal disk is considered. In figure 1, we show the 1^-band 
(5.11 x 10 14 Hz - 6.61 x 10 14 Hz) flux images at different 
binary phases. The mass accretion rates is m = 10 3 , the 
inclination angle is i = 70°, mass ratio is q = Mx/Mq = 1, 
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temperature of the companion star is Tq = 15000 K, and 
velocity of the wind is /3 = 0, respectively. As can be seen 
in figure 1, the geometrical thickness of the disk is thin 
compared with the size of the disk. It is understood that 
the thickness of the disk does not affect the light curve for 
rh = 10 3 . The scale height of the disk weakly depends on 
the rh as can be seen equations (1), (6), (9). 

Figure 2 is the time- variation of F-band flux for vari- 
ous inclination angles. The mass accretion rate is set to 
be rh = 10 3 . When the inclination angle increases, the 
flux drops at phase 0.5. This is because the disk hides 
emission from a secondary star. Optical emission from 
the secondary star is much larger than that from the disk 
for rh = 10 3 . 

As mass accretion rate increases, the disk becomes as 
luminous as the companion star (see figure 3). Some pre- 
vious studies applied a thick disk model to the optical 
light curve analysis in SS433 (Sanbuichi & Fukue 1993; 
Fukue et al. 1997, 1998; Hirai & Fukue 2001). However, 
their model is called "thick torus model" (Abramowicz et 
al. 1978; Madau et al. 1988), and it is thermally /secularly 
unstable (Kato et al. 1998). Apart from the structure of 
the inner most region of the disk, our analytic disk model 
(i.e., slim disk model) is applicable to the supercritically 
accreting regime (Watarai 2006). The model can produce 
more plausible temperature and scale height. Hence, our 
model for a bright object is more realistic than that of 
past researche. 

4-2. Case with Massive Wind 

Figure 4 represents flux images of a model with massive 
wind. The spherical structure at the central region of 
the disk is the photosphere by the massive wind. The 
brightness of the photosphere gradually changes from the 
central region to the outer region. This is called the limb- 
darkling effect. 

Figure 5 shows V-band light curves for high inclina- 
tion angle i = 80°. The most remarkable feature is that 
a primary minimum and a secondary minimum inversion 
happen as the mass accretion rate increases. The F-band 
magnitude at phase or 1 is larger than the magnitude 
at phase 0.5 when the mass accretion rate is relatively 
small (rh = 100 and 1000). However, an inversion of the 
flux happens when rh is large. This is because the optical 
F-band flux from an accretion disk increases as the mass 
accretion rate increases. This feature may be applicable 
to the observation data of SS433. 

In figure 6, we change the wind velocities with various 
values. The size of the wind photosphere depends on the 
velocity of the wind f3. That is, the last scattering sur- 
face is inversely proportional to the velocity (see equation 
(13)). The low f3 outflow therefore makes large photo- 
sphere, and thus the geometrical thickness of the wind 
causes the deep secondary minimum during its eclipse. 
These features are the main results of the present study. 




Fig. 1. V-band flux images in no-outflow models with var- 
ious phases (0, 0.25, 0.5, and 0.75). Horizontal and vertical 
axis are normalized by the binary distance a. The mass ac- 
cretion rate is set to be 1000M cr i t , and mass ratio q = M x /M c 
is q = 1.0. The inclination angle is 70°, and the temperature 
of the companion star T c is 15000 K. 
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0.5 1 1.5 2 

Phase 



Fig. 2. Theoretical V-band light curves expected from our 
disk model (no wind) with different inclination angles. The 
inclination angles are i = 40°,50°, 60°, 70°, and 80° from 
bottom to top. Other parameters are q = Mx/Mc =1.0, 
rh = 10 3 , and T c = 15000K, respectively. 
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Phase 



Fig. 3. Same as figure 2, but as a function of mass ac- 
cretion rates (no wind). The mass accretion rates are 
m = 10 2 ,10 2 5 ,10 3 ,10 3 - 5 , and 10 4 from bottom to top. Other 
parameters are q = 1.0, i = 80° , and T c = 15000K, respectively. 

5. Discussion 

5.1. Eclipsing Light Curves in SS4-33 

One difficulty of the optical light curves in SS433 is 
interpret ating of the secondary minimum at phase 0.5, 
which is made by an eclipse of the companion star by 
the disk. As we explained in the previous section, it is 
difficult to fit the optical light curves observed in SS433 
with only the disk model. The theory also support the 
massive outflow scenario introduced by observations. The 
mass accretion rate in SS433 is appreciably supercritical, 
M ~ 10" 4 M o /yr (van den Heuvel 1980; Shklovskii 1981; 
Perez & Blundell 2009 ). This value is corresponds to 
~ 4.5 x 10 T m for a 1OM black hole, and it seems to be 
an upper limit of the mass accretion rate. 

Since its discovery, binary parameters in SS433 have 
not yet been confirmed. In particular, we could not reach 
consensus as to whether the compact object in SS433 




Fig. 4. V-band images including the effect of wind at phase 
0.5 with various wind velocity (3 = 0.003 (top) and 0.005 (bot- 
tom). Other parameters are q = 1.0, rh = 10 3 , i = 80°, and 
T c = 15000K, respectively. 

is a neutron star or a black hole. Gies et al. (2002) 
evaluated the mass of the companion star via absorp- 
tion lines of A7Ib star to 19±7M . Recently Kubota 
et al. (2009) reported a new constraint on the mass of 
the compact object by the absorption lines taken from 
SUBARU and Gemini observations. They derive the 
mass of the compact object and companion star to be 
M x = 4.ltg'jM and M c = 12.2l§;?M . We apply the 
mass ratio q = Mx/Mq = 0.38 derived by Kubota et al. 
(2009), and the black hole mass M B h set to be 4.OM . 
This mass ratio is not in conflict with other observation 
results (Gies et al. 2002; Hillwig et al. 2004). 

Figure 7 represents V-band images at various binary 
phases. We fix the mass ratio q = 0.38, inclination angle 
i = 78°, and the binary period P = 13.1 day based on 
many former observations. The disk size (radius) is 1.35 x 
10 b r g ~ 1.6 x 10 7 km, which is smaller than the case of 
q = 1.0. The size of the photosphere is comparable to the 
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0.5 1 1.5 2 

Phase 



Fig. 5. V-band light curves with various mass accretion 
rates. The mass accretion rates are ra=100, 1000, 2000, 3000, 
4000, 5000, 6000, 7000, 8000, 9000, and 10000 from bottom 
to top. Other parameters are q = 1.0, (3=0.01, and i = 80°. 




17 0=0.01 



0.5 1 1.5 2 

Phase 

Fig. 6. Same as figure 5, but the dependence of the wind 
velocity (3. The wind velocities are 0.003, 0.005, 0.007, and 
0.01 from top to bottom. 

disk size in the case of ra=5000. 

As for the wind velocity, a quasi-spherical non- 
relativistic wind from accretion disk has 3000 km/s ~ 0.01c 
(Cherepashchuk 2002). Recently, Perez et al. (2009) 
found very fast accretion disk wind by near-IR spec- 
troscopy, and its terminal velocity is about 1500 km s _1 
which is equivalent to 0.5 % of the speed of light. In the 
fit, we fixed the wind velocity to j3 = 0.01 for simplicity, 
allowing changes the other two parameters, namely the 
mass accretion rate and temperature of the donor star. 

Before comparing the observational data in SS433 and 
our model directly, we infer the possible parameter sets 
with figure 8. Figure 8 represents the difference of V- 
magnitude (A My) at phase and phase 0.5 as a function 
of accretion rate. The definition of the A My is given by 

AMy M v (@phase 0.5) - M v (@phase 0). (16) 

In the case of SS433, AMy is about 0.4, which shows 
the magnitude of phase 0.5 is larger than that of phase 
0. Actually, there are some parameter sets that can fit 
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Fig. 7. V-band flux image in SS433 at various phases 
(0, 0.25, 0.5, 0.75). The mass accretion/outflow rate is 
m aC c = m ou t=5000. The black hole mass is Mbh = 4Mq, 
the mass of the companion star is Mq = 12Mq, i.e., mass 
ratio q = Mbh/Mq =0.38. The effective temperature of the 
companion star is Tq =15000 K, and the inclination angle is 
i = 78°, whose values are referred from observational results. 
The velocity of the outflow is fixed at ft = 0.01. 
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Fig. 8. Difference of ^-magnitude in phase and phase 0.5 
as a function of mass accretion rate. Other parameters are 
fixed at M BH = 4M©,M C = 12M©, (3 = 0.01, and % = 78°, 
respectively. 
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Phase 

Fig. 9. V-band light curve of SS433 fitted by our model with 
several parameter sets. Kemp et al. (1986) showed the aver- 
aged optical light curve data of SS433. We set the numerical 
data to agree to the observed flux at phase 0. Other fixed the 
parameters are the same as figure 8. 

the observation, we thus could not uniformly determine 
the best-fit parameter. We tried to fit our model to the 
y-band average light curve by Kemp et al. (1986). It 
is clear in figure 9 that our model can fit whole shape 
of the light curve, but it does not make a distinction 
between three cases, i.e., (i) m=5000, Xc=10000K, (ii) 
ra=7000, T C =12000K, and (hi) ra=10000, T C =14000K, 
respectively. 

The mass-donor (companion) star is supposed to be A- 
type evolved star from the absorption line analysis (Gies 
et al. 2002; Hillwig et al. 2004). On the other hand, some 
recent observations indicate that the companion star has 
lower temperature which is less than 15000 K (Kudritzki 
et al. 2003; Cherepashchuk et al. 2005). It is an impor- 
tant task for observers to confirm the temperature of the 
companion star in SS433. 

SS433 has a precession period of the jet, and the light 
curve changes at each precession phases. It is necessary 



to compare the light curve of each precession phase with 
the model testing. 

In addition, X-ray emission in SS433 comes from not 
only the photosphere of the wind but also the non-thermal 
emission of jet component (Rose 1995; Krivosheyev et al. 
2009; see also Abolmasov et al. 2009). X-ray emission 
from the outflow depends on the acceleration mechanism, 
the temperature distribution, and the initial velocity of 
the outflow launched from the disk surface. If we con- 
sider the X-ray emission more seriously, non-thermal X- 
ray emitting component (maybe jet component) should 
be included in our model (e.g., Reynoso et al. 2008; 
Cherepashchuk et al. 2009). These issues will be our 
future tasks. 

5.2. Comments on Eclipsing ULXs 

Considering the probability of the eclipse events in the 
binary system, it is natural that a few eclipsing binaries 
exist (Pooley & Rappaport 2005). Recently X-ray eclips- 
ing light curves has been detected in several ULXs. As 
the data of the ULXs increases, the number of eclipsing 
ULXs also increases. 

M33 X-7 is one example of an eclipsing black hole X- 
ray binary discovered in recent years, and X-ray eclipse 
has been clearly detected (Pietsch et al. 2004). This ob- 
ject suspects a high mass X-ray binary, but its luminosity 
is very large unlike the famous Cygnus X-l. There are 
several scenarios to explain the high luminosity. Moreno 
Mendez et al. (2008) pointed out the contradiction of the 
black hole spin scenario and proposed a hypercritical ac- 
cretion scenario in M33 X-7 based on the binary evolution 
theory. An eclipsing luminous X-ray binary in the dwarf 
starburst galaxy NGC4214 has been reported by Ghosh 
et al. (2006). They clearly detected the X-ray eclipsing 
feature from this object. 

If both optical and X-ray eclipse are observed in an 
object, we may presume the spatial structure of the ac- 
cretion disk with the difference of the emitting region. We 
are looking forward to waiting for further detection of the 
eclipse in black hole binaries not only in our galaxy, but 
also in galaxies. 

6. Conclusion 

We have calculated the light curves of eclipsing binaries 
by handling more realistic accretion disk models with an 
optically thick outflow. We also applied the present model 
to the supercritical accreting object SS 433. Our calcula- 
tion is somewhat simple, but the geometrical thickness of 
the accretion disk has not been considered seriously so far. 
In addition, we clearly show the change of the shape of the 
light curve by the wind using a numerical calculation. 

As for the model of wind, because it includes unknown 
physics (accretion, geometrical thickness, mass loss rate, 
etc.), so we need to evaluate the temperature at the pho- 
tosphere more seriously. This will be an important issue 
for observed PG quasars (Young et al. 2007). Modeling 
of the detailed physics of outflow is our future issue. 

The measurement of the wind velocity and the mass 
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loss rate via observations of the absorption lines will be 
crucial evidences to confirm our scenario. The fitting with 
the multi-wavelength observational data will be the next 
step to confirming our model. It is also one of our future 
tasks. 
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